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Abstract We present the optical to mid-infrared SEDs of 1 1 debris disk candidates from 
Spitzer SWIRE fields. All these candidates are selected from SWIRE 24/im sources matched 
with both the SDSS star catalog and the 2MASS point source catalog. They show an excess 
in the mid-infrared at 24/im (Ks-[24]vega > 0.44) indicating the presence of a circumstellar 
dust disk. The observed optical spectra show that they are all late type main-sequence stars 
covering the spectral types of FGKM. Their fractional luminosities are well above 5x 10~^, 
even up to the high fractional luminosity of 1 x 10~^. The high galactic latitudes of SWIRE 
fields indicate that most of these candidates could belong to the oldest stars in the thick disk. 
Our results indicate that the high fractional luminosity debris disks could exist in the old 
solar-like star systems, though they are now still quite rare. Their discoveries at high-galactic 
latitudes will also provide us an excellent opportunity to the further studies of properties and 
evolution of the debris disk in the ISM poor environments. 
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1 INTRODUCTION 



One of the most notable achievement of the InfraRed Astronomical Satellite (IRAS) is the discovery of 



dusty circumstellar disks dZuckerman . 



from a main- sequence star (lAumann et al 



200 ih. An example is Vega, which shows a large infrared excess 



19841) . This provided the direct evidence to the existence of a 



debri s disk firstly. The scattered light from several nearby debris disks has been analyzed by use of the 
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disks left over from star formation (IHollenbach et al.l.l200Ql:IWyatt & DentLl20Q2l:lGorlova et al 



study of debris disks is crucial to understand t he formation and ex istence of planets dBryden et al 



smaller objects, such as comets and asteroids (IZuckerman . 



2QQlh . 



2004). The 



20091) and 



Limited by the se nsitivity of IRAS, most of the debr i s disk discovered were for nearby stars much 



younger than the sun (IMetchev et al 



2004 



Brvden et al 



20061) . since they are likely to be at a transi 



tion stage with a h igher frequency to possess a disk with a plenty of dust (IZuckerman & Song . 



Brvden et al 



level ( Krist et al 



2004; 



2006h. an d their luminosities are high enough to heat a debris disk to an IRAS -detectable 



20051). With th e higer sensitivity, better resolution and the more extended wavelength 



(IJourdain de Muizon et al 



searching for circumstellar dust 



(IDecin et al 



2000, 



2003 



1999), the Infrared Space Observatory (ISO) has expanded our capability of 



Spangler et al 



Habing et al 



12001 ) 



20011) to the modest infrared excess among the older stars 



With the launch of the Spitzer Space Telescope (I Werner et al.L |2004|) a new level o f sensitivity and 



spatial resolution is provided to the studies of debris disks in the infrared (ISu et al 



20051) . It provides the 



pote ntial to identify a nd investigate the debris systems that were not detectable wit h previous observato 



ries (Kim et al 



2005k . It also extends the search for debris disks to larger distan ces (IBrvden 



Several 5/? programs have focused on this fie ld, including t he Co res to Disks (lEvans et al 



L et al 



2003 



20061) . 



( Beichman et al 



2005 



Formation and Evolution of Plan etary Systems (IMever et al.L l2004 FEPS), the Galactic Legacy Infrared 



Mid-plane Survey Extr aordinaire ( Benjamin et al.[ 



2003 



C2D), 



GLIMPSE), and the Survey of Solar-Type Stars 



SSS). However, most of these programs are designed for the search of the infrared 
excess in either nearby known planetary systems or the debris systems at regions, where interstellar material 
is rich, such as the low galactic latitude regions and the molecular clouds. While, at high Galac tic latitudes . 



there are also some wide-area surveys with the Multiband I maging Photome t er for 



2004 



MIPS), suc h as the NOAO Deep Wide-Field Survey (IJannuzi & Devi 



Brand et al 



1999; 



Spitzer ([Rieke et al 



Houcketal 



20061) . the Spitzer Wide-area Infrare d Extragalactic Sur vey (ILonsdale et al 



2005 



candidates with infrared excesses (IHovhannisyan et al 



and the Spitzer extragalactic Firs t Look Survey (IFadda et al.L l2006l xFLS), which provide us more star 



2003L SWIRE), 



2009h . 



The Spi t zer SW IRE were observed in the four mid-infrared IRAC bands (3.6, 4.5, 5.8, 8 . fim ) 



( Fazio et al 



2004 and the three mid-to-far infrared MIPS bands (24, 70, 160 /im) 



2004 ). 



Being designed for extragalactic surveys, all the six SWIRE fields are at high galactic latitudes and cover 
about 50 sq. degree. The wide field and high sensitivity prov ide us an opportunity to search for the new faint 
debris disk candidates even the planets dBryden et al.l l2009k in the thick disk or even in the halo. Contrary 
to the other Spitzer programs mentioned above, which focus on the younger or nearby systems, the debris 
disk candidates selected from SWIRE fields are much farther in distance from 300 to 2,000 pc compared 
with those which have a distance from 10 to 150 pc in FEPS, and also much older in age of about 10 Gyr 
belonging to the old est systems. Generall y, the star that is older than 10 Myr can be regarded as a candi- 



date with debris disk dGorlova et al 



20061). This sample will help us to explore the formation of planetary 
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In this paper, we first describe the optical to infrared observations, data reduction and candidate selection 
in Section 2. In Section 3, we give the infrared properties of these 24/im excess stars, their spectral energy 
distributions (SEDs). The discussions and summary are presented in Section 4 and 5. 

2 DATA, OBSERVATIONS AND CANDIDATE SELECTION 
2.1 Spitzer Mid-Infrared Data and Data Reductions 

The 24/im excess stars are from ELAIS-Nl, ELAIS-N2, Lockman Hole fields of the SWIRE fields. The 
common regions of observations between t he IRAC and MIPS for ELAIS-Nl, ELAIS-N2, and Lockman 
Hole fields are about 8,4, and 1 1 sq. degree dSurace et al.Ll2Q05|) . The BCD (Basic Calibrated Data) images 
of the IRAC four bands were obtained from Spitzer S ciences Center (SS C), which include flat-field correc- 



tions, dark subtraction, linearity and flux calibrations (IFazio et al .1120041) . The IRAC images were mosaiced 



from the BCD images after pointi ng refinement , disto rtion correction and c osmic-ray removal with the final 
pixel scale of 0.6'' as described by 



Huang et al 



(l2004 and 



Wu et al 



were mosaiced in the similar way with the fi nal pixel scale of 1 . 225''( Wen et al 



Matching the sources detected by SExtractor (IBertin & Arnouts . 



(120051): Whilst the MIPS 24/im 



2007 



Cao&Wu . 



images 



2007 ). 



19961) in the five bands (IRAC four bands 



and MIPS 24/im band ) with the 2MASS sources, we obtained the astrometric uncertainties of less than 
0.5''. 

The mid-infrared photometries were obtained by SExtractor with an aperture of 3" for t he IRAC four 



bands and 10.0" for the MIPS 24/im band. All these magnitudes are in AB magnitude system dOke & GunnL 



19831). All the magnitudes of four IRAC bands were correcte jd the aperture of 24". Comparing with the 



20031) for the 2MASS stars with 



model colors o f IRAC magnitude s and 2MASS Ks magnitude dCutri et al 

J-Ks < 0.3 as ( Lacy et al]j2005l) did for sources in Spitzer xFLS field, the small additional corrections were 
made for all the four IRAC bands. And these give us the calibration errors in the four IRAC bands better 
than 0.08 mag. As for 24/im, the magnitud es were first correcte d to the aperture of 30". Then an additional 
correction was also performe d according to\l 



Surace et al 



than 10% dRieke et al 



d2005l) . The calibration accuracy of 24/im is better 



20041). The 5a flux in the IRAC four bands and MIPS 24/im band are about 5, 8, 43, 



43 and 200 /iJy respectively. 



2.2 Optical and Near-IR Photometries 

The ELAIS- Nl, Lockman Ho l e and ELAIS-N2 fields are completely or partly covered by the Sloan Digital 
Sky Survey dStoughton et al.l. I2002L SDSS) and the total overlap regions used in this work is about 15 
sq. degree. All the Spitzer sources were match ed with the SDSS point sources from "Star" catalog of 



Data Release 4 ( Adelman-McCarthv et al 



2006h with a radius of 3.0 ". The PSF magnitudes of u,g,r,i,z 



bands were adopte d in this work. The near-infrared magnitudes were from the 2MASS point source catalog 



dCutri et al 



20031). The default J, H, Ks band magnitudes were adopted. 



2.3 Candidate Selection 



Gorlova et al 



( 20041) studied the 24/im emission of the M47 stars, and found that the majority of them 
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they defined the 24/im excess stars as those 3cr redder than the mean value, that is, Ks-[24]vega ^ 0.44. 
We adopted their selection criterion to select 24/im excess stars. [24]ye^a=[24]AB-6.74 was used to trans- 
form the AB magnitudes in the 24/im band to the Vega magnitudes. To obtain the optical spectra, we only 
selected stars with r magnitudes less than 17.5. We also removed the objects with fuzzy features, whose 
24/im emission could come from the background galaxies. Finally, we selected 1 1 24/im excess candidates. 
Figure 5 shows the locations of all the candidates in the color-color diagram and the selection criterion is 
also plotted as a vertical dotted line. Their names and the optical to infrared magnitudes of the candidates 
are listed in Table 1. Since there are problems in the r magnitude of the star Jl 63754. 26+405259.1 and 
the g,r,i magnitudes of the star J163948. 68+41371 1.0, we did not include these measurement in our sub- 
sequence analysis. From Table 1, we can see that except J104508. 69+592830.5, all these candidates have 
Ks-[24]yega valuc abovc 1, indicating strong 24/im excess. The optical r-band, IRAC 3.6/im, 8.0/im and 
MIPS 24/im images are shown in Figure 1. The central circle in each image indicates the optical position 
of each star with a radius of 1.5". It can be seen that all positions in the four bands are consistent. 

2.4 Optical Spectroscopy 

The optical spectra of 1 1 debris disk candidates stars were obtained with the 2.16m telescope at Xinglong, 
NAOC (National Astronomical Observatories, Chinese Academy of Sciences) from February to May, 2006. 
The attached spectrograph is either the OMR spectrograph with dispersion of 200 A/mm or BFOSC (Beijing 
Faint Object Spectrograph and Camera) with the grism G4. Both spectrographs give a resolution of ~ lOA 
and cover the wavelength range from 3 800 A to ~ 8OOOA. The exposure times depend on the apparent 
magnitudes of these stars from 900 to 3600 seconds. The detail informations of the spectral observations are 
presented in Table 2. All these spectral data were reduced by the standard procedures with IRAF packages, 
which include overscan correction for BFOSC only, bias subtraction, flat-field correction. The He/Ne/Ar 
and Fe/Ar lamps were used for the wavelength calibrations of the OMR and BFOSC spectra. The KPNO 
standard stars were obtained to perform the flux calibration at each night. All resulting spectra are shown 
in Figure 2. The spectral classifications are given in Table 2. For those with a poor signal-to-noise spectra 
(Stars 3, 8, 11), we had to determine their spectral types according their SEDs. All these 11 stars are the 
main- sequence dwarf stars with the spectral types from M to F, and belong to solar-like stars. 



2.5 Reddening 



Since all SWIRE regions are at high galactic latitudes, the extincti on by the Milky Way extinction A y of all 



these stars are no larger than 0.04 from the SDSS "stars" catalog (lAdelman-McCarthy et al 



2006) and are 



therefore neglected in the further analysis. How about the local extinction of these stars? Figure 3 shows the 
optical color-color diagram for 9 stars. Due to lack of some optical bands, the stars J163948. 68+41371 1.0 
and Jl 63754. 26 +405259.1 are not plot ted in this figure. As a comparison, the models of main-sequence 
stars and giants (IFitzgerald et al.l.ll970h are also presented as the different curves. All 9 stars are consistent 
with late type stars as above and are located nearby the curve of the model stars. This indicates that all 
these stars are slishtly obscured. The extinction of Kq is only 10% of Ay and that of 24 am is even smaller 
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3 RESULTS AND ANALYSIS 



3.1 IRAC Colors 



Figure 4 shows the mid-IR color-color plot for the 1 1 stars in IRAC bands. The small circles represent the 
positions of the stars and their corresponding numbers are also labelled. Vega is plotted as a star symbol. 
There are not any obvious deviation of colors (either [3.6]-[4.5] or [3.6]-[8.0]) for Stars 1,4,6,9 from Vega. 
As Vega has proved to be an IR excess star, for comparison we also plot the position of black-body with 
Vega temperature of 9600K as plus symbol. Both those stars and Vega show minor excess from black-body 
in [8.0] band. However, Stars 5,7,8 show large deviation (more than 1 magnitude) in color [3.6]-[8.0]. Such 
a 8/im excess can be seen obviously in their SEDs in Figure 6. The remaining Stars 2,3,10,11 show a 
marginal or modest deviation in [3.6]-[8.0], which show small excess at 8/im in SEDs (Figure 6). Only two 
(Stars 3,7) of these stars show a marginal deviation in [3. 6] -[4. 5], and the deviations are about 0.2, about 
2-3 times typical error. Therefore, we can be sure that most of these 24/im excess stars show an excess in 
the 8/im band, and few of them show an obvious excess in the 4.5/im band. 



3.2 MIPS 24/im Excess and Fractional Luminosity 



Figu re 5 presents the J-H versus Ks-[24]vega diagram. The dashed are the curve of main sequence stars 



from 



Gorlovaetal 



(I2004|) . All the 1 1 stars are located at the upper late type stars region. This is consistent 
with the result of the spectral classification. Almost all of these stars present a large 24/im excess. Their 
Ks-[24]vega colors are larger than 1, even to 6. These values are far above our selection criterion value of 
0.44. Only the late type Star 9 has a Ks-[24]vega value of 0.51, little larger than 0.44. 

Unfortunately, we can not obtain the confirmed far-infrared fluxes for all these stars, because of the low 
resolution and low sensitivity of the Spitzer MIPS 70/im and 160/im bands. Without the total infrared lu- 
minosities of stars, we can not obtain the fractional lumin osity fd, which is d efined as the ratio of integrated 



infrared excess of the disk to bolometric luminosity of star 



Moor et al 



dust. To roughly estimate the fractional dust luminosit ies, we assumed Ft r ^ uF^ 



( 2005a ). The bolometric luminosities of stars are from 



(l2006h . to characterize th e amount of 

24/im ], as 



de Jager & Nieuwenhuijzen 



Chen et al 



(119871) . The calculated 

fractional luminosities are listed in Table 2. The fractional luminosities range from 5x 10~^ to 3x 10~^. It 



is quite interesting that four of those stars even have the fractional luminosity values higher than 10 



3.3 Spectral Energy Distributions 



The SEDs of 1 1 stars are shown in Figure 6. They cover the wavelength range from the optical to the 
mid-IR bands, including the available photometries of SDSS u,g,r,i,z, 2MASS J,H,Ks, Spitzer IRAC 
four bands and MIPS 24/im band. These stars exhibit a variety of mid-IR properties. Five of them show 
excess only at 24 fim and four show obviously excess at both 8/im and 24/im. Star 3 even presents excess at 
shorter wavelength of 3.6/im. Except Star 9, all the stars show an devia tion from the photosphere continuum 



at mi d-IR, indicating the existence of an inner hole of the dust disk (IMuzeroUe et al 



20061: 



Young et al 
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3.4 Notes For Individual Source 

J160551.07+534841.0, J160650.59+543420.6, J163754.26+405259.1, J163948.68+413711.0, 
J104205.94+594657.2: These stars have high S/N spectra and a moderate 24/im excess (Fig 5) with 
Ks-[24]vega between 1 and 4. Their SEDs from the optical to IRAC bands are well consistent with 
photosphere continuum and only present infrared excess at 24/im (Fig 6). 

J104508.69+592830.5: This is a M-type star from its spectra (Fig 3). Though it satisfies the Ks-[24]vega 
criterion, it presents a marginal excess above the photospheric emission at 24/im. 

J163236.05+405537.3: This is a G-type star and the only source with 8/im flux higher than 24/im flux in 
the sample (Fig 3 and Fig 6). Such a feature indicates the peak of infrared emission just between 8/im and 
24/im and the outer radius of the dust disk would be much smaller than 100 au. 

J104537. 18+570532.9: Though it has a low S/N spectra, combining with its optical and near-infrared col- 
ors, we still can classify it as a K-type star. Its SED shows that there is little deviation from photospheric 
emission at 8/im. However, it presents a high infrared excess at 24/im (Ks-[24]vega^^)- 
J163611.64+412427.9, J163730.40+403553.1: As above, with the low S/N spectra, J163730.40+403553.1 
was classified as a K-type star with the help of either the colors and the SED. J16361 1.64+412427.9 is a 
G-type star. Both stars show a high infrared excess at 24/im. Their SEDs show an abrupt increase from 
5. 8/im to 8/im and then continue to increase to 24/im. 

J160122.04+545708.2: For this source, we only have very low signal-to-noise spectra and classified it as a 
K-type star. With Ks-[24]vega of 5.62, it is the strongest 24/im excess source among the 11 stars. Its SED 
shows a deviation from photospheric emission even in the near-infrared band, indicating its smaller inner 
radius of dust disk. 

4 DISCUSSION 

4.1 The Coincidence Probability 

Though we excluded the sources with the fuzzy features, it is still possible that some background and 
distant galaxies coincide with the position of the candidates and contribute to the measur ed radiation at 
24/im Therefore, we need to estimate such a coincidence probability for each source. As 



Stauffer et al 



(120051) . we first obtained the cumulative source counts for those with 24/im magnitude brighter than the 
star itself. Here we give an example to explain how we estimated the coincidence probability of the star 
Jl 6055 1.07+53484 1.0. Since it has 24/im magnitude of 16.1, we can estimate the cumulative sources counts 
of 7.8 X 10^ per steradian for sources with the 24/im magnitude less than 16.1 in the SWIRE ELAIS-Nl 
fields based on our catalog. This is consistent with t hat of order of 7x 10 ^ per steradian for flux density 



greater than 1.3 mJy (16.1 magnitude in AB system) (IPapovich et al.Ll2004l) from the another high galactic 



latitude field Spitzer xFLS. It corresponds to about one source per 54,000 arcsec^, considering the previous 
matched radius of 3 arcsec. The coincidence probability of the background 24/im source with magnitude 
less than 16.1 being close to the line of sight to J160551.07+534841.0 is about 0.0005. Similarly, we esti- 
mated the coincidence probabilities of the 1 1 stars are among 0.0003 and 0.003. Since all the stars locate at 
the high galactic latitude, it is impossible to be co ntaminated by cirrus . However, we still can not exclude 
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4.2 Comparison with Other 24/im Excess Samples 



Until now, most detected 24/im excess stars observed by Spitzer are young systems. 



Lowetal 



(2005) 



presented 24/im ob servation for 24 me mbers of the 8-10 Myr old TW Hya association. They found four 



24/im excess stars. 



Young et al 



(120041) detected several 24/im excess stars from young (25 Myr) cluster 



NGC 2547. Almost all of these stars have K s-\24]vpaa colors less than 2. Only one shows a large Ks- 

[24]vega COlor Up tO 4. 



Gorlova et al 



(|2004|) selected 24/im excess stars from 100 Myr old open cluster 
M 47. Seven of the early type stars show smaller excess with Ks-[24]Yega=0.6-0.9 and one early type 
has a color value of 2.4. Three late type stars h ave Ks-[24]Yega colors between 1 to 4. By studying the 



100 Myr Pleiades cluster, 



Gorlova et al 



(20061) obtained an excess fraction of 25% (5/20) for the early- 



type stars and 10% (4/40) for solar-type stars. All of these stars have Ks-[24]vega colors less than 1.5. 



Chen et al 



( 2005ah obtained the observations of 40 F- and G-type members of the Scorpius-Centaurus OB 



association with ages between 5 and 20 Myr at 24/im. They det ected 14 objects that posses 24/im fluxes > 



30% larger than their predicted photosphere. 



Chen et al 



(l2005b|) obtained the observations of 39 A-through 



M-type dwarfs with estimated ages between 12 and 600 Myr. Only three stars possess a 24/im excess. 



Su et al 



(20061) reported the 24/im measurements of 160 A-type main- s equence stars with ages ranging 



from 5 to 850 Myr. The 24/im excess rate is 32%. 



Beichman et al 



(2005) searched for the infrared excess 



debris disks toward 26 FGK field stars known to have one or more plan ets. All these stars have a median 



age of 4 Gyr . None of them show an excesses at 24/im. Including the 



Beichman et al 



Brvden et al 



( 2005l) 's sample. 



(l2006h extended a well-defined sample of 69 FGK ma in-sequence field stars also with the 



median age of 4 Gyr. Only one star shows excess emission at 24/im. 



Beichman et al 



(20061) searched for 



the circumstellar dust around a sample of 88 F-M stars, all detected high S/N 24/im emission, but only a few 



presen t weak 24/im excess, though 12 of them present a significant excess 70 /im emission. 



Morales et al 



(l2005h identified two new d ebris disk candidates with 24/im excess at high galactic latitude SWIRE fields. 



Fajardo-Acosta et al 



latitude xFLS field. 



( 20041) also s earche d for the 24/im excess of main-sequence stars in another high 



Koerner et al 



(120 10|) presented 49 debris disk candidates that were within 25 pc of 
the Sun and with V < 9. The candidate 24/im stars we found in the SWIRE fields are the solar-like stars 
with possible the oldest ages (see the following section). Most of them present large 24/im excess with 
Ks-[24]vega ^vcn approaching 6. Nowadays, only a few main-sequence stars have comparable large 24/im 
excess. 



4.3 The Age 



Though we have obtained the optical spectra of these 24/im excess candidates and classified them based on 
their spectra and SEDs, we c an not determine their ages based on any of the five methods of age estimation 



for the main- sequence stars (ILachaume et al 



19991). The low resolution of our spectra, and even poor S/N 



for some faint stars, prevented us to obtain the ages by either the rotation or the iron abundance methods, 
which need high spectral resolution to measure the rotational velocities or fine absorption lines of different 
metal elements. The absence of calcium emission and the large uncertainty of the distance determination of 
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all these stars are the field stars and not born in a molecular clouds, we can not use the kinematics method, 
too. 

So what we can do is to determine their possible location (disk or halo) in our galaxy. As we know now, 
there are three component s except the bulge in our galaxy. They are the thin di s k, the thick disk and the 



halo. The previ ous works (IBahcall & Soneira . 



1984 



Gilmore . 



1984; 



Du et al 



Ojha et al 



1999; 



Chen et al 



2001 



2003h showed that the scale height of the thin disk were in the range of 24 to 330 pc and the 



Du et al 



(120061) analyzed 21 BATC 



scale height of the thick disk is in the range of 580 to 1300 pc. Recently,!] 
fields with 15 intermediate-band filters. The scale height they obtained for the thin disk varies from 220 to 
320 pc and those for the thick disk varies from 600 to 1100 pc, which is consistent with previous results. 
Considering that all the three SWIRE fields are at high galactic latitudes of about 40 degree to 50 degree, 
more than half of the candidate stars have a vertical distance to galactic plane above 600 pc (estimated from 
their spectral types), which comparable to the scale height of the thick disk. Therefore, these stars could 
belong to the thick disk or even the halo. The thick disk is an old component with an age around 10 Gyr 



( Gilmore et al 



ISSslllSSsi. It is quite possible that many of these candidates belong to the population of 



the very old stars. 



4.4 The Fractional Luminosity 

There are two major zones of debris in the solar system: the asteroid belt at 2—4 au composed of rocky 
material that is ground up by collisions to produce most of the zodiacal dust clo ud and the Kuipe r Belt 



(KB) that consists of small bodies orbiting beyond Neptune's orbit at 30 — 50 au (IKim et al 



2005). The 



excess emission at the mid-IR band Spitzer 24 /am and IRAS 25 /im are sensitive to material at several au 



dSu et al 



20061). Therefore, the stars with 24 /im excess are probably those with an pla netary systems 



containing a planetesimal belt corresponding to the aste roidal zone of the so l ar sys tem (IGorlova et al 



Backman & Parescd (Il993|) and 



Dermott et al 



20061) . From combination of observation and modelling, 
(20021) determined that the fractional luminosity of our asteroid belt is 10~^ to 10"''. 

Artymowiczl (Il996h argued that debris disks are confined to fd < 10~^ and sources with higher 



fractional luminosity 
"transition" objects). 



probably contain a 



signifi cant amount of gas (e.g. T Tau and Herbig Ae/Be stars. 



Zuckerman & SongI (120041) hypothesized that the stars with fd>lO ^ are younger 



than 100 Myr, and therefore a high fd value can als o be used as an age indicator. This is supported by 



Observations with the ISO by 



Beichman et al 



(|2005|) . who suggested a decline in the fraction of stars with 
the excess IR emission with time, but still with the possibilit y of existing modest (fd > 10~^) excesses 



among the older stars (IDecin et al 



2000 



20031) . Meanwhile, 



Brvden et al 



(2006) selected stars without 



regard to their age, metallicity, or any previous detection of IR excess. The stars have a median age of ^4 
Gyr. Their result shows that the debris disks with > 10~^ are rare around old FGK stars. 

On the contrary. 



Decinetal 



(|2003|) claim ed the existence of high fd disks around older stars. Based 



on the general evolutionary trend described by 



Moor et al 



(|2006|) . the fractional lumin osities of individua l 
systems were found to show a large spread ofl0~^ < fd < 10~^at almost any age (IDecin et al 



20031) . 



Particularly interesting is the relatively high number of older systems (t> 500 Myr) with the high values of 
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Most of our candidates have the values of the fractional luminosity between 10~^ to 10~^ and four 
of them even present the high fractional luminosities greater than 10 All the four stars have a vertical 
distance to the disk above 900 pc and have the high possibility to be the old solar-like stars in either the 
thick disk or in the halo. Therefore, they are possibly the rare systems up to date. 



4.5 Physical Explanation and Challenge 

How such a large amount of dust still exists in very old solar-like stars? Due to the effect of radia- 
tion pressur e, Poynting-Robertson drag and collisional destruction, the lifetime of the dust gra i ns are 



quite short (IBeichman et al 



Dominik & Decinl. 



2005 



Rieke et al 



2005 



Backman & Parescd. 



1993 



Lagrange et al 



200C; 



20031), for example, the dust grain with the smaller size (< 1/im) has a blown out time 



less than 100 yrs and t he larger dust gra in drifted by Poynting-Robertson drag would be destroyed on the 
timescale of 1-10 Myr (IKim et al.Ll2005h . These are all far shorter than the a g e of the central star s. Therefore 



the dust observed must have been recently produced (IBeichman et al.U2005l) . lRieke et al.l (l2005h pointed out 



that these reproduced dust in the "debris disk" would arise primarily from collisions between planetesimals 
and from cometary activity. However, the current model of gas planet formation and isotopic evidence from 
terrestrial and lunar samples in dicate that Jupiter mass planets and Earth-Moon sy stem would formed within 



several tens millions of years (ISilverstone et al 



2006; 



Kleine et al 



2002 . 



20031) . This still can not explain 



the high fractional luminosity ph enomenon of the o ld stars. They raise a serious challenge even to the new 



generation of theoretical models (IMoor et al 



20061). The most likely explanation for the presence of debris 



disk with the high fractional luminosity at a ges well above Gyr is dela y onset of collisional cascades by late 



planet formation further away from the star (IDominik & Decin . 



20031). However whether such a mechanism 



can also explain the very old star with an age of 10 Gyr is still questionable. 

Besides, the most important is to confirm the ages of these stars with the high fractional luminosity. 
Therefore the higher S/N and higher resolution spectra are needed in the future observations. 



5 SUMMARY 

We present 1 1 debris disk candidates with distance from 300 to 2000pc and age about lOGyr from Spitzer 
SWIRE fields, with the help of SDSS star catalog and 2MASS point source catalog. They all present 24/im 
excess. The optical spectra obtained from 2.16 telescope at Xinglong, NAOC and SEDs from the optical to 
mid-infrared are also presented. All these observations shows that: 

All of these stars present the late type spectra and the infrared fractional luminosities from 5 x 10~^ to 
3x 10"^. They would be solar-like stars with the planetary debris disks. Except J104508.69+592830.5, the 
SEDs of all other stars indicate the existence of an inner hole of the dust disk. 

We infer that many of these candidates would have ages of 10 Gyrs, since they are located in the thick 
disk or in the halo of our galaxy. Therefore, four of these stars could belong to the oldest stars with the 
high fractional luminosities (/^ > 1 x 10~^). Though they are quite rare until now, we indicate that the high 
fractional luminosity debris disks can exist in the old solar-like star systems. 

With the release of Wide-field Infrared Survey Explorer (WISE) data, there will be more and more 
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latitudes will provide us an opportunity to study the properties and evolution of the debris disk evolution in 
ISM poor environments. 
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Fig. 1 The SDSS r-band, Spitzer 3.6, 8.0 and 24/im bands images of the 24/im excess stars. 
The circle in each image gives the position of stars in r-band with a radius of 1.5". 
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Fig. 1 Continued. 
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Fig. 2 The observed low resolution optical spectra of the 24/im excess stars. Star 3, 8 and 11 
have very low S/N. All the stars present the late type star features. 
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Fig. 3 The optical B-V versus U-B diagram of the 24/im excess stars. Each star is symbolized 
as small circle and labelled with the number. Stars 4 and 6 are lack, because of problems in 
some optical bands. The dotted line shows the normal dwarf stars in the color-color plot, and the 
corresponding spectral types are also labelled. 
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Fig. 4 The mid-infrared IRAC color-color diagram of the 24/im excess stars. All the colors are in 
AB magnitude system. Each star is symbolized as small circle and labelled with the number. The 
star and plus symbols present the positions of Vega and black-body with temperature of 9600K. 
Three stars present a high mid-infrared excess at 8/im band. Eour present a modest excess and 
the other four show a minor excess at 8/im band. 



18 



Wu H. et al. 



X 

I 




Ks-[24] 



Vega 



Fig. 5 The plot of J-H versus Ks-[24]vega- All the color are in the Vega magnitude system. 
Each star is symbolized as small circle and labelled with the number. The dashed curve shows 
the normal dwarf stars labelled with the corresponding spectral types. The dotted line gives our 
criteria to select the 24/im excess source. Except Star 9, all the stars have Ks-[24]vega higher 
than 1 magnitude and four stars present very strong 24/im excess. 
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Fig. 6 The optical to mid-infrared SEDs of the 24/im excess stars. The black points give the 
fluxes of the SDSS five bands (u,g,r,i,z), the 2MASS three bands (J,H,Ks) and the Spitzer 
IRAC four bands and MIPS 24/im band. The photospheric emission is ploted as solid curves. 
Comparing with the solid curves, all these stars present excess at mid-infrared. 



Table 1 The names, magnitudes and colors of the 1 1 24/im excess stars. 



NO. 


Name 


u 


9 


r 


i 


z 


J 


H 


Ks 


[3.6]as 


[4.5]as 


[5.8]as 


[8.0]as 


[24]as 


Ks-[2A]vega 


1 


J160551.07+534841.0 


14.64 


13.13 


12.52 


13.38 


12.29 


11.41 


11.02 


10.95 


13.71 


14.25 


14.65 


15.29 


16.13 


1.56 






±0.02 


±0.00 


±0.00 


±0.00 


±0.01 


±0.02 


±0.02 


±0.02 


±0.08 


±0.08 


±0.08 


±0.08 


±0.10 


±0.10 


2 


J160650.59+543420.6 


19.39 


16.83 


15.52 


14.88 


14.56 


13.35 


12.67 


12.52 


15.13 


15.61 


16.07 


16.42 


15.94 


3.32 






±0.04 


±0.02 


±0.02 


±0.02 


±0.02 


±0.03 


±0.02 


±0.02 


±0.08 


±0.08 


±0.08 


±0.08 


±0.10 


±0.10 


3 


J160122.04+545708.2 


18.84 


17.46 


16.86 


16.69 


16.66 


15.78 


15.38 


15.12 


17.68 


17.96 


18.12 


18.56 


16.24 


5.62 






±0.03 


±0.03 


±0.02 


±0.02 


±0.02 


±0.07 


±0.11 


±0.13 


±0.08 


±0.08 


±0.08 


±0.09 


±0.10 


±0.17 


4 


J163754.26+405259.1 


16.10 


14.94 


- 


14.39 


14.37 


13.49 


13.16 


13.10 


15.86 


16.26 


16.72 


17.40 


16.86 


2.98 






±0.02 


±0.01 


±0.00 


±0.01 


±0.02 


±0.02 


±0.03 


±0.03 


±0.08 


±0.08 


±0.08 


±0.08 


±0.10 


±0.11 


5 


J163236.05+405537.3 


16.34 


15.28 


14.85 


14.76 


14.77 


13.82 


13.54 


13.43 


16.13 


16.58 


16.96 


16.39 


16.41 


3.76 






±0.02 


±0.01 


±0.01 


±0.02 


±0.02 


±0.03 


±0.05 


±0.05 


±0.08 


±0.08 


±0.08 


±0.08 


±0.10 


±0.11 


6 


J163948.68+413711.0 


14.85 


- 


- 


- 


13.25 


12.44 


12.19 


12.18 


14.90 


15.35 


15.88 


16.51 


17.01 


1.91 






±0.01 


±0.05 


±0.03 


±0.03 


±0.01 


±0.02 


±0.02 


±0.02 


±0.08 


±0.08 


±0.08 


±0.08 


±0.10 


±0.11 


7 


J16361 1.64+412427.9 


17.81 


16.37 


15.82 


15.67 


15.60 


14.73 


14.27 


14.18 


17.02 


17.27 


17.72 


17.14 


15.83 


5.09 






±0.02 


±0.02 


±0.01 


±0.02 


±0.02 


±0.03 


±0.05 


±0.06 


±0.08 


±0.08 


±0.08 


±0.08 


±0.10 


±0.12 


8 


J163730.40+403553.1 


19.38 


17.94 


17.28 


17.06 


16.93 


15.90 


15.30 


15.36 


17.74 


18.06 


18.46 


18.01 


16.97 


5.13 






±0.03 


±0.01 


±0.01 


±0.01 


±0.02 


±0.08 


±0.10 


±0.19 


±0.08 


±0.08 


±0.09 


±0.08 


±0.10 


±0.22 


9 


J104508.69+592830.5 


18.88 


16.08 


14.69 


13.60 


13.02 


11.77 


11.09 


10.89 


13.45 


13.83 


14.32 


14.99 


17.11 


0.51 






±0.05 


±0.01 


±0.02 


±0.02 


±0.01 


±0.03 


±0.03 


±0.02 


±0.08 


±0.08 


±0.08 


±0.08 


±0.11 


±0.11 


10 


J104205.94+594657.2 


16.61 


15.45 


15.08 


14.97 


14.93 


14.11 


13.83 


13.70 


16.57 


16.92 


17.46 


17.83 


17.23 


3.21 






±0.01 


±0.03 


±0.02 


±0.01 


±0.02 


±0.03 


±0.04 


±0.04 


±0.08 


±0.08 


±0.08 


±0.08 


±0.11 


±0.12 


11 


J104537.18+570532.9 


20.30 


18.12 


17.21 


16.90 


16.69 


15.64 


15.04 


15.02 


17.53 


17.88 


18.34 


18.60 


16.75 


5.01 






±0.05 


±0.02 


±0.02 


±0.02 


±0.02 


±0.07 


±0.11 


±0.14 


±0.08 


±0.08 


±0.09 


±0.09 


±0.10 


±0.18 
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Table 2 The log of observation, spectral type and the fractional luminosity of the 11 24/im 
excess stars. 



NO. 


Name 


Date of Obs. 


Exptime 


Instrument 


Slit 


Sp. 


fd 








sec 




arcsec 






1 


J160551.07+534841.0 


Feb. 4, 2006 


900 


OMR 200A/mm 


2.5 


G8V 


7.0e-5 


2 


J160650.59+543420.6 


Apr. 28, 2006 


2400 


BFOSC Grism#4 


1.8 


MOV 


7.4e-4 


3 


J160122.04+545708.2 


May 6, 2006 


3600 


OMR 200A/mm 


2.0 


KOV 


3.2e-3 


4 


J163754.26+405259.1 


May 5, 2006 


2700 


OMR 200A/mm 


2.0 


G2V 


2.1e-4 


5 


J163236.05+405537.3 


Feb. 21, 2006 


3600 


BFOSC Grism#4 


1.1 


G3V 


4.2e-4 


6 


J163948.68+413711.0 


Feb. 4, 2006 


1800 


OMR 200A/mm 


2.5 


F8V 


5.4e-5 


7 


J16361 1.64+412427.9 


May 5, 2006 


2700 


OMR 200A/mm 


2.0 


G7V 


1.8e-3 


8 


J163730.40+403553.1 


May 6, 2006 


3600 


OMR 200A/mm 


2.0 


K4V 


2.1e-3 


9 


J104508.69+592830.5 


Feb. 3, 2006 


2400 


OMR 200A/mm 


2.5 


M3V 


6.9e-5 


10 


J104205.94+594657.2 


Feb. 3, 2006 


2400 


OMR 200A/mm 


2.5 


GOV 


2.4e-4 


11 


J104537.18+570532.9 


Feb. 3, 2006 


3600 


OMR 200 A/mm 


2.5 


K6V 


2.3e-3 



